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Abstract 

We study a quite general family of dynamical r-matrices for an auxiliary loop algebra C(su(2)) 
related to restricted flows for equations of the KdV type. This underlying r-matrix structure 
allows to reconstruct Lax representations and to find variables of separation for a wide set of the 
integrable natural Hamiltonian systems. As an example, we discuss the Henon-Heiles system and 
a quartic system of two degrees of freedom in detail. 
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1 Introduction 



The main aim of this paper is to study the r-matrix formulation for mechanical systems embedded as 
restricted flows into KdV type equations and to investigate its connection with Lax representations 
and with the method of separation of variables. Lax representations are essentual for solving these 
mechanical systems either by linearization on the Jacobian of determinant curve of the Lax matrix 
Jul or by determining variables of separation through the functional Bethe ansatz |27], |2{| . 

Restricted flows are easiest understood as stationary flows of soliton hierarchies with sources ||. 
Many well known integrable mechanical systems, such as Henon-Heiles, Gamier, Neumann (and many 
other) can be embedded into the KdV and other soliton hierarchies as restricted flows j6| 

The main advantage of such embedding is that important structure elements for restricted flows, 
such as Lax representation, bi-Hamiltonian formulation, Newton parametrization |^4| can be system- 
atically derived from the underlying soliton hierarchies. But it is not the case, as yet, for the r-matrix 
formulation of these mechanical systems. There is no direct way of transporting the well known r- 
matrix formulation for soliton hierarchies |p^| to their stationary and restricted flows. The underlying 
r-matrix structure for finite-dimensional systems appears to be more complicated then for integrable 
PDE's. A new type of, so called, dynamical r-matrices depending both on the spectral parameter and 
on dynamical variables has to be considered. The first example of such r-matrix has been found for 
the Calogero-Moser system j^o) and then for the parabolic and elliptic separable potentials jy], [?2| . 

It appears that the parabolic and elliptic r-matrices belong to a quite general family of dynam- 
ical r-matrices which we introduce in this paper through a convenient ansatz which generalizes our 
r-matrix from jl^, [52|. We study the general algebraic properties of these r-matrices and derive Lax 
representations through proper specializations. These specializations are guided by the known inte- 
grable cases of natural Hamiltonian systems related to restricted flows for equations of the KdV type 
|^|, ^8). Further we use these Lax representations for finding variables of separation an explain how 
starting with variables of separation these Lax representations be reconstructed by using the Sklyanin 
approach [§7|, |§. 

The paper is organized as follows. 

In section 2 we consider additive algebraic deformations of linear r-matrix algebras on an auxiliary 
loop algebra g — C(su(2)). Kernels of these r-matrices have either rational, trigonometric and elliptic 
dependence on the spectral parameter [|l2| [l3[ ^5|. We consider deformations L(X) = Lq(X) + F(X) 6 
£(su(2)) of a certain Lax matrix Lq(X) € £(su(2)), which satisfies the standard r-matrix brackets. 
These deformed Lax matrices obey a dynamical rs-matrix algebra. From dynamical point of view the 
deformed Lax matrices correspond to systems with additively shifted integrals of motion 

rnew jold , r 

4 = 4 + tk- 

Further for r-matrices of the XXX and XX Z type we consider projection of deformations F(X) on 
finite-dimensional Poisson (adjj-invariant) subspaces Cm,n(su(2)). Their Lax matrices Lmn{X) obey 
the rs-matrix brackets as well. At the end of section 2 we construct Lax pairs from the rs-matrix 
algebra and prove that these Lax pairs are related to stationary flows for equations of the KdV type. 

In section 3 we express matrix elements of the matrices L (X) and Lmn(X), defined on a direct 
sum of loop algebras ®£(su(2)), through canonically conjugated variables of separation and prove 
that these matrices obey the rs-matrix brackets. Such matrices describe a geodesic motion and a 
potential motion on a Riemannian manifold with complex diagonal meromorphic Riemannian metrics 
of special type. Corresponding potentials are the finite gap potentials on these manifolds. The problem 
of quantization is discussed briefly. 

In section 4 we consider finite-dimensional integrable systems of natural type and investigate three 
types of canonical transformations of variables prescribed by the linear r-matrix structure. Two of 
these are pure coordinate transformations to generalized elliptic coordinates and to Jacobi coordinates, 
while the third one provides an example of a momentum dependent change of variables. We exemplify 
of our approach with the Henon-Heiles system and a quartic potential of two degrees of freedom. 

In the concluding remarks we discuss shortly an application of our constructions to other nonlinear 
integrable PDE's such as the AKNS and the sine-Gordon hierarchies. 
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2 Deformations of linear r-matrix algebra 



2.1 Additive deformations 



Let us recall basic algebraic constructions in the r-matrix scheme 25, |26| for loop algebras, which we 
shall use in this paper. A loop algebra g = C(a) is defined as follows. Let a be a Lie algebra, then 
C(a) consists of Laurent polynomials X(X) — J2 x j^ with coefficients in a for which commutator is 
defined as [x\ t ,y\ : >] = {x,y]X l+: >. There is an obvious decomposition of g = £(a) into a linear sum of 
two subalgebras 

g = g++g-, g+=®j>oaA 3 , g-=($ j<0 a\ 3 . (2.1) 

Let P± be projection operator onto g± parallel to g T . A standard i?-matrix associated with £(a) is 
defined as 

R = P+ - P_ . (2.2) 
We shall use tensor form of i?-brackets which is more common in the inverse scattering method p2l 



|13| , p.7[ . If the algebra a is identified with its dual then the standard R- matrix (2.2) is the Hilbert 
transform in C(a) with the kernel 

r(A,/i) = T ^- ) (2.3) 

where t 6 a ® a is a kernel of the identity operator in a. The adjoi nt k ernel is r*(A, /i) = — r(A, /i). For 
twisted subalgebras Q{a) e £(a) the corresponding decomposition ( |2.l| ) yields kernels with trigonomet- 
ric and elliptic dependence on the spectral parameter jl2| |l3|, pjjj|| . In the quantum inverse scattering 
method r-matrices with rational, trigonometric and elliptic dependence on the spectral parameter are 



called r-matrices of XXX, XX Z and XYZ types respectively |12|, |13|, [17 . 

In this paper we consider the Lie algebra a = su(2) only. We fix also representation of the algebra 
su(2) and shall work in the matrix notation, although the underlying constructions are independent 
of particular matrix representation. In the following we shall consider spin one half representations of 
su(2) or shall use the direct sum of these representations. This means that we restrict ourselves to 
two-dimensional auxiliary space in which su(2) matrices act. 

There are two important Lie-Poisson brackets (i?-brackets) |6[ ||l[] which can be consedered 
as linear classical limit of the fundamental commutator relations in the quantum inverse scattering 
method ||, |§ . The y 

arc the r-bracket 



{L(A),1( M )} = [r(A, M ),L(A) +1( M ) ] , (2.4) 



and the rs-bracket 



{L(X)Mp)} = [r(A, (J,),L(\) +L(p) } + [s(X, /i),l(A) -L(ji) } , (2.5) 

l 

which is related to the reflection equations pq, B2|. Here we use the standard notations L( A) = 

2 ^ 

L(X) ® /, L([i) =7® L{p). The brackets ([2.4|^.5]) are Lie-Poisson brackets, if the matrices r(A, ^l) 
and s(A, fj,) satisfy the modified classical Yang-Baxter equations (2^, |2^, Q. 
For a traceless 2x2 matrix Lq (A) (which we shall call Lax matrix) 

3 

L (A) =5> fc (A)-a fc (2.6) 

k=l 



which obeys the linear r-matrix algebra (2.4) the r-matrix is an antisymmetric function of one argument 



r(A, fi) = r(A — /i) and r(A) = — r(— A). Namely, put 

(2.7) 



k=l 



The Ok are Pauli matrices and coefficients Wk{X) are functions of spectral parameter only. They are 
determined by construction of R- matrix on loop algebra C{a ) (pJ |2~2 ) p5| , or by the requirement 
that the r-matrix obeys classical Yang-Baxter equations 
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It follows from the algebra fl2.4|) that the determinant e?o(A) = det Lq(X) can be taken as a generating 
function of commuting integrals of motion I k since 

{d (X),d ( fJ ,)} = 0, d (A)=]T/ fe A\ A,M6C (2.8) 

fc 

implies {Ik, Ij} = . 



A similar property is valid for the rs-algebra ( |2.5| ) too. 

Here we consider deformations L(X) £ £(a) of the Lax matrix L {X) £ £(a) with a = sm(2). The 
main stimulus for considering such deformations has been inspired by the study of Lax representation 
and r-matrix formulation for restricted flows of the KdV and coupled KdV hierarchies ||, Q, which 
describe many physically interesting integrable natural Hamiltonian systems such as: Neumann, Gar- 
nier, Henon-Heiles systems and an infinite family of integrable polynomial potentials |2j| . The relevant 
deformations of Lq(X) had the form 

L{\) = ( n(h ° \) (X)=L (X) + F(b,X) (2.9) 



G(b,X) + c -a 

where matrix L Q (X) £ £(a) and hence the entries a(A) , 6(A) and c(A) are Laurent series satisfying 
Poisson brackets ( |2.4| ) with the r-matrix of the XXX type. The deformation G(6, A) had a special 
form connected to KdV equation [jl8[ [52| . 

Here we find a natural generalization of this result where the deformed matrix L(X) belongs the 
same phase space. The additive deformation of £o(A) has the form 

= £ (s k (X) + a*(A) ■ 4^-) • o k = L Q (X) + F(X) , (2.10) 
k=i V 26(A)/ 

3 

where 6(A) = £ a»(A)s 4 (A) , (2.11) 

i=l 

and a\(X) + a 2 2 {X) + a 2 3 (X) = 0. (2.12) 

Here functions aj(X) and the ar bitrar y function /(A) are complex- valued functions of the spectral 
parameter A only. The condition ( 2.12| ) guarantes that determinant of L(X) has the additive property 

d(X) = d (A) + /(A) , 

which means that integrals of the motion d(X) = IkX k of the deformed integrable systems differ by 
some constants 

h = h + Jk, Jk e ^ • 



The requirement that the deformed Lax matrix L(X) satisfy the rs-bracket (2.5) with the same 
r-matrix r(A — //) yields four algebraic equations for the unknown functions afc(A) , k = 1, 2, 3 and a 
certain function g(X,fi). 

Theorem 1 The deformed Lax matrix ( 2. 1 satisfies the linear rs-algebra (2.L ), if 

w j (X,^)a j (ii)a i (X) - w l {X,n)a i {^)a j (X) = g{X, fi) ■ a k (X) , (2-13) 

where (j, i, k) are cyclic permutations of indices (1, 2, 3) and the scale function g(X, /i) depends only on 
the spectral parameters X and fi. The corresponding matrix s(X, /i) is given by 

3 

s(X,fi)= £ aij(X,[i) oi ® a 3 , (2-14) 



with 



a, (A, fi) = Wt (X,») • 9 (X,») • _ WjM . g{/h A) . ^(X)f(X) 

b z (n) o (A) 
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Proof: The proof is a direct but lengthy computation. 

Notice that the function /(A) is quite arbitrary here. A compact form of the equation (2.13) follows 
from the ideas of Sklyanin in p7[ : 

Lemma 1 The condition ( 2.1^ ) of Theorem^ for the coefficients ccfc(A) is equivalent to equation 

{6(A), Kjjl)} = g(X, fx) ■ 6(A) - g(fx, A) • 6( M ) . (2.16) 



Proof: This lemma is proved by direct substitution of the equality (2.13) into the equation (2.16) and 
vice versa. 

Let us present some of the best known jl2[ |l7) r-matrices and the related solutions for ctk(X) 
and g(X, (x). The r-matrix of the XXX and XX Z types are 



Wl = W2 = w(X) 



n 



9 {X) 



and W3(\) 



with 
and 



A 

sinh A 



in the XXX case , 
in the XX Z case , 



(2.17) 
(2.18) 



here rj is constant and ip'(X) denotes derivative with respect to A. 

An interesting particular class of solutions to equations (2.12- 2~13| ) is distinguished by the condi- 
tions 



a 3 = . 



ai = . 



(2.19) 



They represent two points on the circle in complex plane with the phase difference equal to tt/2, for 
instance ai = 1 , «2 = i- For these solutions the function g(X, fi) is equal to zero {g(X, fi) — 0). A 
different example of solutions for r-matrix of XX Z type is 



ai = cosh A , a2 — i sinh A . 

Usually one uses the following rational parametrization 
for r-matrix of XXZ type: 



013 



w\{u — v) = U)2(u — v) = rj- 



2uv 



W3(u — v) = Tf- 



(2.20) 



1 

u 



with 



as, for an example, in the description of the sine-Gordon equation [ fL2[ 
An elliptic r-matrix of the XYZ type has 



with 



u>i(A) = r\ 

W 3 (A) = T) 

a% = -1 ; 



9 / 11 9 10 (A,fc) 
ei e n (A,fc) ' 

Qji e 01 (A,fc) 
6oi e n (A,fc) ' 

i dn(A, k) 
k cn(A, k) 



w 2 (X) = r) 



e , 11 9 o(A,fc) 
eooe n (A,fc) 



(2.21) 



tt2 = 



a 3 



A; cn(A, k) 



Here fc and fc' denote modulus of the Jacobi elliptic functions and 0^ -(A, k) are the elliptic theta function 
in the notation of p4j. We have to emphasize that the functions 6(A) (2.11) with the coefficients afe(A) 



(2.21) and the function g(A,/x) were introduced in [|27j for the quadratic r-matrix algebra as solution 
of the equation (2.16). The function g(A,/i) is independent from second parameter /x function and it 
is equal to 

9j x e w (X-K,k) 
Bio 0ii (A — K, k) 
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We recall from p7] that, for the purpose of separation variables, Sklyanin considered a matrix Lq(\) 
which is similar to the Lax matrix Lq(A): 



with 



L (X) = 



U(X) 



a b 

c —a 



©oi @oo 



(X) = U- 1 (X)L (X)U(X) : 
2£ 



(2.22) 



1 + fc 



2y/k 
1 + k 



Here matrix Lo(X) obeys the linear r-matrix algebra ( |2.4| ) with the r-matrix ( [2.21 ). The similar matrix 
Lo(X) satisfies the r-bracket (2.4) with r-matrix of the form 



p(A,/x) = U- 1 (X)U- 1 (X)r(X,f,)U(X)U(X) 

3 

= 22 w j (A — /i, A;) • (Tj (g) <Tj + 

3=1 

+ [w 3 (fi- K,k) -w 3 (X- K, k)} ■<r 3 (g)a 3 + 
+ iw(p — K, k) ■ a 3 <g) <J2 — iw(X — K,k) ■ o~i 



' cr 3 ; 



(2.23) 



w(X — [i, k) 
w 3 (X - fi) 



= V: 



wi(X - fi, k) — w 2 (X - /1, k) = 77 

e' n (A- M ,fc) 



e , 11 e 10 (A- / i,fc) 

O10 6n(A - fj,, k) 



6n(A - fi, k) ' 

here if is a standard elliptic integral of the module k. Notice that now two first coefficients of the 
matrix p(X) are equal wi(X) = W2(X) in same way as it was in the XXX and XX Z cases. 

It is well known that r-brackets (^J) is the Lie-Poisson brackets if r-matrix obeys the classical 
Yang-Baxter equations []l2|, ||(| . The Yang-Baxter equation is an equation on matrices in the space 
C 2 (£) C 2 <g) C 2 . We use the familiar tensor notations Lj for the matrix acting as matrix L in the j's 
factor of the product and trivially acting in the remaining factors, for example L\ = L ® I ® I . In 
similar way a matrix r%j is acting trivially in the third A:'s factor and it works as an matrix r in the 
product of the other r/'s factors. Classical Yang-Baxter equation for pure numerical r-matrices has 
the form 

[d±(A, (i), df 3 (X, 1/)] + [df 2 (X, (i), 4 3 ( M , u)] + [4,(1/, /i), ^3 (A, 1/)] = . (2.24) 

Following [llj we introduce dynamical Yang-Baxter equations for the dynamical rs-matrix. Let us 
consider matrices = r ± s corresponding to particular class of solutions ( [2.19 ). 

Theorem 2 For the rs-algebra (2.1 ) with matrices (2.7- 2.ll\ ) the following equations are valid 



[df 2 {X,fi),df 3 (X,v)} + [d± 2 {X,fi),d± 3 (fi,v)} + [df 2 (v, fj,), df 3 (X, v)] + 
+[L 2 (ix),df 3 (X,u)} - [L 3 (i/),4(A,M)] + 



+[X(X,(,,v),L 2 (n)-L 3 (v)} = 0. 
The matrix X(X, /x, v) is 

X{X,\i,v) — (3(X, [i, v) ■ O- <£> er_ ® cr_ , 



(2.25) 



(2.26) 



P(X,p,u) 



f(X)b- 3 (x) + /(m)6- 3 (m) • iy - A) + ■ (A - M ) 



(A - ^t)(/i - ^)(i/- A) 

T/ie other two equations are obtained from ( 2.2£\ ) by cyclic permutation. 

Proof: This Theorem also can be proved by a straightforward computation. 

Tensor Y(A,/i,i/) is completely symmetric tensor with respect to any permutations of auxiliary 
spaces. For the quite arbitrary solutions au to the equations ( |2.12^T3| ) we can introduce an asymmet- 
ric tensors X^ tk ' (A, fj,, v). Then dynamical Yang-Baxter equations take the general form introduced 
in ft 



G 



2.2 Special deformations of the XXX and XX Z cases and Lax representa- 
tion 

This section is based on list of formulae collected in jy], [lj, |32| and on meaning its in general r-matrix 
approach to loop algebras. 

Our purpose is to apply the rs-matrix formulation for description of natural Hamiltonian systems 
in order to solve them through separation of variables. For separating variables we shall apply the 
functional Bethe ansatz ^7], . The Bethe ansatz has an asymmetric matrix form (l^, [2!| (variables 
of separation are introduced as zeroes of the entry 6(A)). Therefore we restrict ourselves to a particular 
class of deformations obtained by certain projections of the matrix 



L(X) = 



f(X)b~ l +c 



(A) = L (A) + 





/(A)fo- 1 



(2.27) 



where the entries a(A) , 6(A) and c(A) are absolutely convergent Laurent series in the XXX case or 
absolutely convergent Fourier series in the XX Z case 



a W = J2k a k exp(/c ■ A) 



for XXX case, 
for XX Z case, 



and similarly for 6(A) and c(A). 

It is well known [£5| that the standard -R-bracket on £(a)* associated with ( [2.l| [^2]) has a large 
collection of finite-dimensional Poisson (ad^j-invariant) subspaces 

C m ,n = ®f=- M a * xj , provided M > ; N > -1 . (2.28) 

In other words the subspaces Cm,o & n d Ci,n are invariant under the coajoint action of the subalgebras 
g + and g_ (|2~l|). 

Let us introduce the following projection operators [-]m7V on these Poisson subspaces 



\z\mn 



\z\mn 



J2 ^x k 



N 



E **** 



MN 



k=-M 



(2.29) 



+ OC 



^ z k exp(fc • A) 



k=-c 



N 



z k exp(k ■ A) . 



J MN 



-M 



for the Laurent and for the Fourier series. In particular in the XXX case 



£ z k X k 



+ 00 



]>>A fc , 



(2.30) 



k=0 



denotes the Taylor projection. 

We shall define a new additive deformation of L (X) as 



Lmn(X) = L (X) + F MN (b, A) 



a b 

Ga/at(6, A) + c -a 



(A), 



where Ga/./v(6, A) is a function of the spectral parameter A and of the entry 6(A). It reads 

G M N(b,X) = [f(X)-b- 1 (X)] MN , 

oo oo 



k— — oc 

For the Taylor projection this deformation takes the form 



k= — oo 



N 



G N (b,X)^[f N (X)-b- 1 (X)}. 



where /jv(A) =^fk A fe 



(2.31) 

(2.32) 
(2.33) 

(2.34) 



fc=0 
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The essential feature of this deformation is, in comparison with ( 2.10 ), that the determinant of the 
modified matrix Lmn(X) is different from do(X) + f(X) 

d(A) = do (A) + b ■ G MN ^ d (A) + /(A) 



and therefore integrals of motion I% ew of the deformed system are functionally different from the 



undeformed integrals I£ ld , and yet the Lmn(X) matrix belongs to a certain rs-algebra. 



Theorem 3 The matrix Lmw(A) ( 2.31 ) satisfies the linear rs-matrix algebra (2J.), with the matrix 
saw(A,/i) given by 

s(\,n) =a MN (A,/i)cr_ ®<7_, er_ = ( !j* ^ 
TTie junction olmn (A, /i) is defined by 

a MN (A,/*) = -KA-M) ([/(A)^ 2 (A)] M jv - [/(/z)&~ 2 (m)]m;v) , (2.35) 
where w(X) = u>i(A) = 102(A). 

Proof: The proof is a straightforward calculation. 

The Theorem 3 has a transparent meaning in the framework of the general r-matrix approach to 
l oop algebras. The Theorem 3 asserts that the matrix s(A, /1) can be restricted to subspaces £mm 
( [2.28 ) and subspaces £m,n is a common Poisson subspace for the rs-brackets (2.5) with matrix sain = 



[s(A, /i)] m n ■ Since the deformations (|2.6| ) are directly connected with the KdV equation |18j, Theorem 
3 describes the relationship between compatible Poisson structures of the KdV equation and the 
Virasoro algebra [pf. Matrices d ± — r ± smn obey dynamical Yang-Baxter equations, which are 



obtained by restriction of the equation (2.25- 2.26|) to the corresponding subspace Cmn- For certain 



matrices Lmtv(A) and smn this proposition has been proved in [ JLjJ . 

In order to describe Lax representations for integrable Hamiltonian systems related to matrices 
L (X) and Lmn{X) we introduce matrices [^5[ 

di2(A,/i) — r(X,fi) + s(A, /i) , and <2 21 (A,/i) = r(X,(J.) — s(A, /i). 

In terms of these matrices the Lie-Poisson rs-brackets reads: 

{£(A),I( M )} - [di2(A,/i), L(A)] + [d 2 i(X,fi), L( M )] , (2.36) 
For constructing Lax representations we shall make use of the following lemma ||] . 



Lemma 2 If matrix L(X) obeys the Lie-Poisson algebra (2.3t ) then the spectral invariants 

H n (X)=tv(L n ), 7i = I,2,...;, (2.37) 
are integrals of motion in involution {H ni H m } = and the Lax equation are given by 

L( M ) = {H n (X),L(»)} = [M n (A,/i),L(/i)], (2.38) 



1 



with M n (X,fi) =ntri ^""^A) d 2 i(A, /i) J , n= 1,2... . (2.39) 

where tri means trace in the first auxiliary space and dot over L means derivative with respect to time 
corresponding to the Hamiltonian H n (X). 

In two dimensional auxiliary space 

d(X) = detL(X) = -i(tr (L 2 ) = ~H 2 (X) . 

In order to get a Hamiltonian H, which does not depend on the spectral parameter, one needs a 
linear functional (a projection) which selects, for example (in XXX case), a coefficient at certain 
power of A 

H=l$ x [H 2 (X)]=$ x [d(X)]. (2.40) 
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For instance it can be defined as a residue of order m at Ao 

1 ef™" 1 



$ A [z] = ResZ z(X) = 



(2.41) 



A=A 



Then the Lax representation for the Lax matrices (2.6 or 2.27) with the Hamiltonian ( |2.40 ) has the 
form 



Up) 



{H, L(p)} = {-$ A [tn L 2 (X)],L(fi)} = [M(jj), L(jj)] 



(2.42) 



M{ji) = 



tn L(A)-dai(A,/*) 



= $ A [m(a,m): 



In particular for the initial Lax matrix Lq (2.6) one gets s — 0; c?2i(A, /i) = r(A — fi) and then 

L (/z) = [M (n),L Q (p)} (2.43) 



fc=i 



M (Ai) = $ A 
In this paper we start with the requirement that 

M = o+ = 



2^ s fc (A) • w fe (A - ^t) • cr fe 



1 




(2.44) 



(2.45) 



and choose L (X) and $a appropriately. It is a rather strong restriction for the matrices Lq (2.6) and 
for the corresponding Hamiltonians. This requirement comes from the study of rs-representations for 
natural Hamiltonian systems related to stationary flows for equations of KdV type |l8|, |32]|. 

As an immediate consequence of the Lax representation ( 2.42] ) with the matrix Mo ( 2.45| ) we can 
rewrite the Lax matrix Lq in the form 



? b 

J xx u x 

T ~2 



b x = {H Q ,b}, 



(2.46) 



where Hq is a Hamiltonian corresponding to Lq- Since the determinant of Lo(A) is a generating function 
of integrals of motion we see that it obeys the equation 



{H , d Q (\)} = -d , x (X) =dl-b = b xxx = . 



(2.47) 



Now for a fixed triad {Lq, Mo, $a) we introduce a modified matrix L(X) (2.27). In this case the 
matrix M(/i) is constructed according to ( 2.42| ) with c?2i (A, /i) = r — s and with the use of the same 
projector $\. It reads 



M( M ) = 



1 

-u(jjl) 



u(M) = f(»)b- 2 (fi) . 



(2.48) 



due to the linearity of <!>a 

$a [«* ■ ( -/(A) 6- 1 (A) + c(A) + (/(A)6- 2 (A) - f i^ 2 {»)) ■ b(X) )] = 

= -f( f ,)b- 2 ((i) = ~u( f i). 



(2.49) 



Then from the Lax representation (2.42) it follows that the deformed matrix L(X) has to be equal to 



L(X) 



-b(X)u(X) 



, b x = {H,b}. 



(2.50) 
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where the Hamiltonian H — $\[d(X)] is obtained by applying the old projector <3>a to the modified 
determinant. The determinant d(X) of the modified matrix L(X) obeys the equation 



-d x (X) 



ud x + -u x ) ■ b = Bi[u] ■ b = 0. 



(2.51) 



Here B\ [u] is the Hamiltonian pencil oper ator for the coupled KdV equation J?], ||] . 

For the deformed (according to (2.31)) matrix Lmn{^X) the Lax equation is constructed analo- 
gously (2.48-2.50) and (2.51). One has to replace u(/i) = /(/i)& _2 (/i) with the potential Umn(h) = 



\MN- 



Equation (2.51) shows that our deformation is connected to the stationary problems for equation 
of the KdV type. The Hamiltonian and some properties of motion for the finite-dimensional integrable 
systems connected with stationary flows of the KdV equations were considered in || ^, ||. Lax 
representations, bi-Hamiltonian structures, Newton representations and some other properties for these 
systems were introducing in J?], |^, 24 through method of the restricted flows for nonlinear equations. 
Some particular systems were studied in the r-matrix approach JTl| [l8[ |32] | . In the next section we tie 
these different methods to the Lie algebraic method of constructing matrices L(X) on the loop algebras. 

In the next section we shall assume that the Lax matrix Mq has the form 



Mn = 



*'() 





Uq 
-Vq 



(2.52) 



It introduces certain technical complications without influence on the conceptual structure. 



3 Construction of Lax matrices L(X) for separable systems 
3.1 Separation of variables 

The main problem for a given integrable Hamiltonian system with a complete set of functionally inde- 
pendent and commuting integrals of motion Hi , . . . , H n is to determine its solutions. The construction 
in the Liouville theorem has a local character and there is no general way to describe solutions globally 
unless some additional information about the system is known, such as a complete, spectral parameter 
dependent Lax representation or variables of separation. 



Here we study the Lax matrices satisfying the linear rs-algebra (2.4-2.5) and method of separation 
of variables. By separation of variables we mean here construction of n pairs of canonically conjugate 
variables 

{v,j,Uk} = {vj,v k } = 0, {vj,Uk} = Sjk , j,k = l,...,n (3.1) 
and functions such that 

*k(uk,v k ,H 1 ,H 2 ,...,H n )=Q, k = l,...,n (3.2) 

where Hk are integrals of motion in involution fl29| . Note that the canonical transformation from 
the original variables (xfe, pk), k = 1, . . . , n to variables of separation (iik, i>k), k = 1, . . . , n may not 
necessarily be a pure coordinate change. It can involve both coordinates and momenta. 

For the Lax matrices L(X) satisfying r-matrix algebra with the r-matrix of the XXX type the 
problem of determining variables of separation has been essentially solved by Sklyanin (see Theorem 
4 below). 

Here we provide a solution of the converse problem how to construct Lax matrices L(X) satisfying 
the rs-matrix algebra for an integrable system given in terms of variables of separation. We express 
first L(X) in terms of variables of separation (v,k, Vk), k = 1, . . . ,n and then we shall use the canon- 
ical transformations naturally determining from the r-matrix structure to the "physical" variables 
(x fc , p k ), k = l,...,n. 

Theorem 4 (Sklyanin [ ^7| , p9[ |) Coordinates Vk, k = 1,. .. ,n defined as 

b(X = u k ) =0 v k — a(X = M fe ) (3.3) 
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are canonically conjugate variables of separation, if the matrix L(X) obeys the r-matrix algebra with 
r -matrices (2.11) or (2.18). Separation equations are 



v l = -d(u k ) , 



(3.4) 



where d(A) = detL(A). 



Proof: We repeat the Sklyanin proof for the XXX and XXZ cases, because this technique will be 
used later. 

Let us list commutation relations for b(fj.) and a(A), 



(3-5) 
(3.6) 



{6(A), b(n)} = {a(A),o(/*)} = 0, 
{a(A), &(/*)} = w 3 (X-n)b^)-w(X-^)b(\), 

where we denote that w(\ — /i) = ti?i(A — (i) = W2 (A — /i). The e qua lities {it;, Uj} — follow from ( |3.5| ). 
To derive the equality {ui,Vj} = we substitute u = Uj in (3.6), and obtain 

{Vj,b(ji)} = w 3 (uj - n) b(ji) . 

It together with the equation 

= {vj,b(ui)} = {vj,b(fj,)} | M = Ui +b'(ui){vj,Ui}, 

gives 

{vj,Ui} = _ ^-y{ u i' 6 (^)} U=«*= 

Equalities Uj} = can be verified in a similar way: 

-{v i: Vj} = {a(ui),a(uj)} 

= {a(uj),a(fx)} | M=Uj +a / (ui){«i,o(ttj)} 
= a'(uj){a(uj),Ui)} + a'(ui){ui, a(y,j)} = 0. 

After the substituting A = Uk into the identity 

a 2 (A) + 6(A)c(A) = -d(A) 

we obtain separation equations 

u| = -d(«fc), fc=l, ...,n., 

where d(A) is the determinant of L(X). 

In the XFZ case we have to use entries a(A) and 6(A) of the similar matrix L(X) ( 2.22 ). By the 
Lemma |l| th e zeroes of entry 6(A) are commuting variables {ttj, uj} = 0. For the rest of the proof we 
refer to p7fl. 



3.2 Construction of Lax matrices L(A) 

A general construction of matrices Lo(A) with multiple poles on a direct sum of loop algebras © C (su(2)) 
has been presented in detail in p5[ . The Lax matrix is assumed in the form 

3 N 

io(A) = ^^s> a (A-^)<7 a , (3.7) 

a=l j 

where <7 a are Pauli matrices, w a (X) are coefficients of the corresponding r-matrix and Vj are simple 
poles of the fixed divisor D — {v\, • • ■ , ^jv}- Residues s" are the standard linear coordinates on su(2)* 
with Lie-Poisson brackets 



11 



. The corresponding phase space is the direct sum of N copies su(2)*. These matrices describe sys tem s 
of N interacting Euler tops. The case of higher order of poles can be treated as a degeneration of ( |3.7| ). 
It describes different tops as well. 

We shall express matrix elements of Lq(X) through canonically conjugate variables. Let us consider 
mcromorphic matrix- function L (X) and the coefficients u>i(A — fi) = W2(A — fi) of the r-matrices 
(2.17-2.18) with one distinguished point at infinity (it is zero by both arguments and their residues at 
infinity are equal to —1). Let us choose the linear functional <&>, as the higher residue fixed order K 
at infinity 



$x[z] = -ReO(A) . 



Then (2.42 and 2.52) gives 



-Res 



K 



6(A) 



= u ■ 



-Res 



K 



-Res 



K 



if{X - fj,) 

c(A) 
<p(X - n) 
q(A)^(A-/i) 

<p(X - n) 



0. 



v ■ 



(3.8) 

(3.9) 
(3.10) 
(3.11) 



In particular uq = 1 , vq = for the matrix Mq ( 2.45 ). In the XXX case it means that entries 6(A) and 
a(A) have the highest in A terms in power K + 1 equal to u and vq, respectively. The e ntry c(A) has 
the highest term of or der K only. Note that the conditions for the entry a(A) and c(A) ( 3.1 0| - 3.11 ) in 
the Lax pair (Lq, Mq) ( 2.45 ) follow from the first condition (3.S) for entry 6(A) and from the definitions 
a(A) = —b x /2 = —d x b/2 and c(A) = —b xx /2, since the derivative d x commutes with the functional $\. 

According to the construction in J25| we can fix one divisor D consisting of poles in the matrix- 
function Lq(X) and a second divisor U consisting of zeroes of the entry 6(A). Let they have the form 



D = {(ej,mj) j = l,...,m;}, 
U = {(uk,n k ) k = l,...,n;} . 
u fe ^ oo , and ej ^ u k , Vj, k . 



(3.12) 



where D, U C CPi = CU {oo} for the XXX model and D, U C CPi/2ttZ for the XX Z model. We 
consider the simple zeroes u kl n k = 1; k = 1, . . . ,n, only. This restriction is related to the r-matrix 
structure. The poles at ej can be multiple poles and their orders rrij are fixed. In the formulae below 
each pole is counted according to its multiplicity and we don't use special indices when it is clear from 
a context. 

Now we introduce the following ansatz for the entry 6(A) according to (2^, [27j 



6(A) 

<p(\) 
<p(X) 



u 



w-^X) = X 
= sinh A 



(3.13) 



in the XXX case , 
in the XX Z case , 



where uo is either a constant or a dynamical variables. The function ip(X) depends from the type of 
the r-matrix and the functional &\ ( |3.8| ) is residue with K = 1 for n < m and K = n — m + 1 for 
n > m. 

The entry a(A) has m poles at points ej of the divisor T> and its values at n points u k are given 
(a(A = Uk) = Vk) ( [3.3| ). The asymptotic behavior at infinity is fixed by Lax representation ( |3.9| ). By 
Lagrange interpolation formula the entry a(A) can be represented as 



a(X) = 6(A) • lv + ^ Vi ■ w 3 (X - u t ) ■ gi, 



where 





db 




9ii = 


dX 


A — u i 



= Res Mi 6 _1 (A) = 



(3.14) 
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The w 3 (X) = ip'(X)/(p(X) is the coefficient in the definition of the corresponding r-matrix. Variables 
Mo j are connected to the Lax matrix Mq (2.52). 

Further, from the Lax representation (2.46) for the initial matrix Lq(X), we can construct the entry 
o(A). It reads as 



o(A) 



— 7T- ■ 1" > Ui ■ W 3 (X - Ui) 

2 V u ° ti / 



where 



iik = Uk,x = d x Uk , and 103(A) 



(3.15) 
(3.16) 



We set now, f or th e simplicity, uq = 1 and vo = as in the Lax matrix ( 2.45| ). Definition of entry c(A) 
follows from ( 2.46 ) as well 



/ \ \ "xx 



a 2 (A) 6(A) 



6(A) 



2 ' ^ 

»=1 



w 3 {X - Ui) —u t 7^- 



(3.17) 



where iij and Ui are defined below. 

From comparison of two forms of the entry a)X) (3.14) and (3.15) we can introduce the following 
Hamiltonian 



»=1 i=l 

using the Hamilton- Jacobi equations. 

The remaining symbols in the entry c(A) are defined as 

Ui = — — = 2 • i>i • gu , 



"96(A) 




<9A 





(3.18) 



■ 9u ■ 



3=1 



if{ui - Bj) 



+2 ^2 [i v k + Vi) 2 gugkk - vf{g kk - gu)gu] ^(u l - u k )(p'(ui - u k ) , 



fe=i 



where y(A) = u) x X (A) = A or <p(X) = sinh(A). Then a direct calculation shows that the Hamiltonian 
( [3.18D is equal to the Hamiltonian @.4C| ) H = $ A [ ri (A) ]. 

The matrix Lo(A) with entries ( |3 . 1 3| — |3.14 ) and (3.17) obeys the r-matrix structure (2.4) as we can 
check it directly. Let us calculate 



{a(A), &(/*)} = w 3 (X-(i)b(n)-w(X-n)b(X) = 



2^ , /1 rx-, t ,.w„-,„.r 6 ( A ^ = ^L^' 



-J y>(A - u k )ip(p - u k ) ^ \cp(X-u k ) ip(ii-Uk)J 

-1 re— 1 

_ y> gfcfc • - ^fc) y> l'(A - ApyQi - - y(A - Uk) 

^ ip{X - u k )(p{fi - Uk) ¥>(X - n)<p(\ - Uk)y{n - Uk) 

The last equlity is true due to the additive theorem for <p(X). 

So, we have constructed the initial matrix Lq(X) in the term of variables of separation. From the 



initial matrices Lq(X) we can construct modified matrices Lmn(X) (2.31). These matrices describe 
integrable Hamiltonian system with the Hamiltonian H = $>\[d(X)] and a complete set of the integrals 
of motion can be obtained from determinants of these matrices as well. The determinant c?(A) is a 
meromorphic function with the divisors of poles determined by divisors of the entries of Lq(X) or of 
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Lmn(X). Therefore the rational algebraic functions d(X) can be written as a quotient of two product 
or it can be decomposed into simple fraction according to the Mittag-Lefner theorem 



d(X) = 



Q7W n;=i ^(a - e 3 ) 



(3.19) 



N—m 



fc=l 



where N = n for the initial matrix Lq(X). Both decomposition provide us with sets of integrals of 
motion. 

Lemma 3 Zeroes hk of the function Pj^(A) can be chosen as a first set of integrals of motion in 
involution for the Hamiltonian systems related to the Lax matrices Lq(X) and Lmn{^)- Zeroes 
of the function Pj v_m (A) and the residues Hji at points Cj constitute a second set of the integrals of 
motion in involution. 



Proof: The proof of involutivity is based on the equality fl2.8p and follows the same pattern as the 
proof of Theorem |j. Completeness and functional independence of these integrals have been proved 
by Mishchenko and Fomenko (see review p5[). 
We summarize these considerations as 



Proposition 1 The Lax matrix Lq(X) (3. IS, 3.14, 3.11) describes geodesic motions on the Riemannian 
manifold with the complex diagonal meromorphic Riemannian metrics 



i kj = si ■ 



"96(A) 




dX 


X=Ui 



-x> n^Li ^ - e i) 



(3.20) 



defined on moduli of n- dimensional Jacobi varieties in terms of complex elliptic (root) coordinates u^. 
The Hamiltonian system (3. IS) 



has a complete system of first integrals Vk and it defines a Lagrangian submanifold of the phase space 
C 2n which has the form of symmetric product ((T x . . . x T)/a n ) of n copies of Riemannian surfaces 



& 



d(A)|, 



Q4X) 



(3.21) 



\=u k 



A modified Lax matrix Lmn(X) (2.31) describes potential motion on the same Riemannian manifold 
with a Hamiltonian of natural type 



H = H Q + V M n (ui,u 2 ...,u n ; f- M , f- 



M+l, 



. In) ■ 



(3.22) 



fj are coefficients of the function /(A) = fkX k or of /(A) = X) These matrices obey the 

rs-matrix structure (2.1 ) with the smn -matrix given by the Theorem^. 

In the potential case we reconstruct, for instance, motions under permutationals symmetric potential 
classically (via point transformation) separable in generalized elliptic coordinates ]23] . For m < n + 1 
we have here the well known Jacobi problem of geodesic motions on quadrics |3| |7|, corresponding 
quadrics are defined by condition 



-Res 



K 



6(A) 



ip(X - \i) 



= 1 



with variables Xj will be introduced in section 4.1. 
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For these systems we can define a Riemann matrix of periods of holomorphic differentials and by 
using the Abel-Jacobi map we can find the action-angle variables. The potentials Vmn are the finite 
gap potentials on these manifolds jlO|. |) . For details about geodesic and potential motion with metric 
( ^20|) in the XXX case see §, |, §f 

Thus we have constructed initial matrices L (X) and modified matrices Lmn 

(A) ([^-JOI]) which 

correspond to a geodesic motion and a potential motion on Riemannian manifolds with metri c (|3.20 ), 
respectively. For the XXX case the initial matrix L (X) ( j3 . 13| , [3.14 . [3.17 ) was considered in |16| and 
the Hamiltonian Hq ( 3.1 8| ) describes free motion on Riemannian spaces of constant curvature. 

The above results describe motions on Jacobi variety ((r x . . . x T)/a n ) which is a symmetric 
product of uniform Riemannian surfaces 



r 



d(A)| ; 



(3.23) 



This scheme can be generalized to the case of the more general Jacobi variety ((ri x Tj . . . x r„)) with 
different Riemannian surfaces 



<i) . „(i)2 



<p)(x) 



1. 



1, 



\=Uk 



,n. 



(3.24) 



where fc and j are two numbers indexing curves. Motion on such manifolds can be described with the 
use of the block matrices 

L(X) = 0™ L«(A, v{ j \u{»; . . . ; t,g>, , (3.25) 

acting in the auxiliary space 

V aw = ® j V<g>,. (3.26) 
Here j numbers blocks and k denotes number of degree of freedom related to each block. Each block 



LW)(A,«P,uW;...;t,g),uW) 



is a matrix of type L (X) or LmnW ( P-27| ), respectively. These manifolds have been investigated in 
p3| and such block matrices were considered in |l5|, ^] . 

In the XYZ case parameterization of the entry 6(A) is more complicated. For instance, it has to 
comply with the equality ( 2.16 ) 



{6(A),6( M )}= ff (A,M)-6(A)- ff (/x,A)-6( M ). 

In order to bypass this difficulty we can use two-dimensional lattice averaging ]l2]| for the constructing 
matrices L (X) from the corresponding matrices of the rational model. 

In quantum mechanics the metric (3.20) defines standard operator of the Laplace-Beltrami type 
and stationary Schrodinger equations 



V J 'Vj* + (E- V MN (u u u 2 ..., u n ; /_ 



M ■ 



Jn))* = 



on the n-dimensional Riemannian manifolds ||] . Quantum separation equation for geodesic motion 
and potential motion and their properties have been considered in |4j || but not in the r-matrix 
approach. In the framework of the linear r-matrix formalism separation of variables has a direct 
quantum counterpart which is described in 1 27, 28, E9|. In the quantum case we replace variables 



Ui, Vi by operators and Poisson brackets ( |2.4| ) by the commutators 



L(A),L( M ) 



= [r(A,/x),L(A) +L(m)! 



Note that dynamical smn(X, /i) matrices depend on coordinates only. Separated equations have the 
following form 



+ U{u k ; e\,... ,e m ; f_ M , . . . , f N ; hx, ■ ■ ■ , h n ) ) ipki^k, hx, ■ ■ -,h n ) = 
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where potential U(uj) is a rational algebraic function of the variables Uj parametrized by a, i = 1 ... to; 
and fj, j = —M . . . N and hj are the eigenvalues integrals of motion. According to psfl , common 
multiparticle eigenfunction of quantum integrals of motion is equal to 

n 

*(it x , . . . ,u n ) = Y[ ipk (wfe) ■ 

k=l 

A quantum operator Lq(X) for free motion in the XXX case was investigated in [^8[ [l9| while, the 
quantum modified operator L(A) for potential motion was considered in fll| , [l8| . 

4 Canonical transformations of variables 

We can represent meromorphic matrix-functions L(X) by means of generators of various Lie groups, 
which one can consider as the various phase spaces for the integrable systems. This leads to links among 
different integrable systems: the Neumann problem, Gaudin magnets and Euler tops p5| , |l6)| . In this 
paper we shall consider finite-dimensional integrable systems of natural type and shall investigate three 



types of canonical transformations of variables prescribed by linear r-matrix structure (2.4.2 



4.1 Generalized elliptic coordinates 

If the limit orders of poles rrij at the points ej are fixed then the Poisson brackets at different poles 
are independent and the corresponding orbits are a direct product of orbits of algebra su(2). This 
factorization is indepen dent of the poi nts e 7 |25| . By the Mittag-Leffler theorem we can decompose 
rational functions 6(A) (|3 .13 ) and <z(A) ( 3.14 ) into partial fractions 



6(A) = pn-m (A) + bji ■ - ej ) = [] - A) + £ hi ■ <^(A - e 3 ) 



k=l 



(4.1) 



o(A) - Pr m - 1 W+Y^a ji -v- i (\-e j )= ]J <p{a k - A) + ^ aji ■ ^(A - ej ) , 

j,i k=l ji 



where bk and are zeroes of the polynomials P"~ m (A) and P" _m_1 (A), respectively. The residues 
bji and aji at the points ej are linear coordinates on the algebra su{2)* (3.7). We use later two 
representations for generators of this algebra in terms of canonical variables 

b 3 = x) , aj = XjPj , Cj = v) , (4.2) 

bj = 2xjXj , a,j = XjPj + XJ'j , Cj = 2 Pj Pj . (4.3) 

We omit here the double indexes ji which are necessary for higher order poles and we used notation 
Cj for residues of entry c(A). Variables Xj , Pj and Xj ,Pj are canonically conjugate coordinates and 



momenta. Generalized elliptic coordinates correspond to the first representation (4.2) and in our 
scheme they are a simple consequence of the Mittag-Leffler theorem for meromorphic functions. 

The other (n — m) variables can be chosen as zeroes of the polynomial P fc "~ m (A = Xk) =0 and as 
the values of function P" _m_1 (A) at these zeroes Xk 

P b n - m (\ = x k ) = p k = P^- m - 1 (X = x k ) (4.4) 

Variables x k , k = m + 1, . . . ,n are calculated from the coefficients of the polynomial P^ _m_ (A) = 
J2BkX k . For instance in the XXX case with simple poles these coefficients are equal to 

.2 



fc-i 



p k+i 



then the conjugate (to the x k ) momenta pk, k = to + 1, . . . , n are 

Pk 



3 = 1 J 
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Lemma 4 Variables 
conjugate variables. 



,Pj J 



= 1, 



and Xk,Pk k = m + 1, . ..,n are a system of canonically 



Proof The Mittag-Lcffler theorem defines a unique, invertible map between old n pairs of variables 
(itfc, Vk), k — 1, . . . , n and the new n pairs of variables (xk,Pk), k — 1, . . . , n. The proof of canonical 
conjugation of these variables is based on the r-matrix structure. The technique developed by Sklyanin 
(see |27| and proof of the Theorem ||) has to be employed. 

In the case of higher order poles we have an additional freedom in representing residues through 
canonical variables. There are two approaches to this problem. The first one is connected with standard 
degenerations of simple divisors, for example ej — > or ej — > oo. This approach is developed in 
fief [l9| . Our approach is different. 

We shall show how the variables are connected to Newton type variables for stationary 

flows of KdV type equations introduced in |24| . We shall consider the XXX case but generalization 
for the XX Z case is transparent. The r-matrix algebra is associative; a linear combination of any two 
matrices L\ and L2 satisfying (^4) satisfies (2.4) too. Let us consider the p olynomial part of Lq(X) 
corresponding to poles at infinity. The determinants of Lq(X) and L(X) from ([2.47-2.51) are 



do(X) 



b-K 



4 ' 



(4.7) 



d(X) 



b ■ b x 



- -2 +b 2 -umn(X) 



They are integrated forms of the KdV recursion relations j24j . 
A simple substitution for the entries of matrix Lq(X) 

6(A) = B 2 , a(X) = -b x /2 = -BB x . 
c(A) = -b xx /2 = -Bl-BB xx , 



(4.8) 



(4.9) 



turns determinants (4.7-4.8) into the form 

rfo(A) = B 3 B XX , d(X) = B 3 B XX + B 4 



/(A) 



B 1 



(4.10) 



MN 



if we use an explicit formula for the potential umn(X). These equations have the form of Newton 
equations for B 



B x ^ 



B,. 



d(X)B 



-3 



B 



d a (X)B- 3 . 
/(A)- 



(4.11) 



B 4 



MN 



If we assume that B = X^^Lo Q N -j A J is a polynomial then its coefficients qj obey the Newton equation 



of motion ( 4.11 ) with d(X) — X/fcA fc , where Ik are integrals of motion. It becomes clear in terms of 
variables qj that deformation of integrals of motion Ik affects only the potential (g-dependent) part. 
The kinetic (momentum dependent) part of Ik remains unchanged. 

Coefficients of all entries a(X) , 6(A) and c(A) are easily expressed in terms of coefficients of B. Let 
us fix the highest order of A in entry as K, then 



K 



K 



6(A) b, X? 



t (A) = ]Ta^ 



K 

c(X) = J2^ j , 



with bx = 1 , ax = , ck = 

due to the definition of M and <I>a- Let 



(4.12) 



(4.13) 



K 



S (A) = ^2 1k-jX j , with q = 1 
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where the set B we consider as a formal set. After substitution of these definitions into ( |4.9| ) we obtain 

K-j K-j 



i=0 



IK- 



E 

i=0 



Qi.x IK- 



(4.14) 



E 9i > x 



i=0 



E 

i=0 



: QK-j-i 



where q x is derivate with r espec t to x and we use d the Newton formulae for a product of sets. From 
qo = 1 and the definitions ( 4.14 ) the restrictions ( |4.13 ) for the higher coefficients of entries of matrix 
Lq follow automatically. Canonically conjugate (to the coordinate s q.j) momenta pj can be derived 
from the r- matrix algebra (2.4) and from the definitions ( 4.12 - 4.14 ) 



Pj — QK+i-j,x i 
As an example, we present first polynomials 



{Pj,Qk} = 8. 



jk 



K 
K 



K = 3. 



6(A) 
6(A) 

-o(A) 
6(A) 

-a(A) 



A + 2gi , 
A 2 + 2Axi 4 
Ap2 + (pi + 
A 3 + 2A 2 xi 
A 2 p 3 + A(p 2 



-a(A) =pi , 
(2x 2 + a; 2 ), 

P2X\) , , 

f A(2x 2 +x 2 )- 
+ P3^l) + (Pi 



(4.15) 



2(x 3 - 
-P2X1 



xix 2 ) , 
-P3X2) ■ 



All integrals of motion can be expressed in terms of new variables directly from the equations ( 4.1C ) 
by taking residues. Notice that the kinetic part of the Hamiltonian H — $^[d(A)] has the nondiagonal 
form 

A'+l 

(4.16) 



T = V PjPK-j 



3=1 



Definition of variables of separation is not changed by deforming the Lax matrix Lq(A) ( p. 27 ). For 
instance, the free motion Hamiltonian H = T (4.16) and the corresponding Hamiltonian with potential 
terms H = T + Vmn separate in the same coordinate systems. Later we shall show some examples 
illustrating the above scheme. 



4.2 Jacobi coordinates 

Let us take a simplest matrix Lq (A) with the entry 6(A) of ( |3 . 13 ) which has only one zero 



L fe) (X, v k ,u k ) = 



-v k <p(X - u k ) 
v k 



It obeys the linear r-matrix algebra ( [Q| ) with the r-matrix ( ^.17| ) or ( 2.18 ). Let 

L (A) =®lL<£\\,v k ,u k ) 



(4.17) 



(4.18) 



what means that we associate to our system a matrix L, which has blocks structure and each block 

(k) 

L obeys the r-algebra with the same r-matrix. 
A modified matrix L(X) (2.27) is equal to 

L(\) = ® k l L$ N (\,v k ,u k ), 

where blocks L^ N (X,v k ,u k ) arc 



(4.19) 



L MN(^ v k, u k) 



- V k 



/(A) 



ip(X - ufe) 



MN 




(4.20) 



—Vk f{X - u fe ) 

CMN(X,U k ) V k 
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These matrices are deformations of the initial matrices Lo(A, v k ,Uk) (4.17) and each block L^ N obeys 
the rs-algebra with a common r-matrix but different (since they depend on Uk) matrices s^ N con- 
structed according to ( 2.35 ). 

The matrix Lq (4.17) and the modified matrix Lmn (4.20) have a internal structure pl| 



L (X,v k ,u k ) 



LMN{\v k ,Uk) = 



-Vk <p(X - u k ) 
Vk 



(fc) 



r (fc) 

E OLjP) 



- E a jP 

/(A) 



(fc) 



¥>(A-E«7 1 ?j fc) ; 



^(A-E^^f) 



E«iP 



(A:) 



MN 



if variables and itfc are considered as linear combinations of certain canonical variables pj^ , j = 
1,. . .K. If we consider a system with of n degree of freedom and require that Hamiltonian has a 
natural canonical form with the kinetic energy form 



(4.21) 



k=l 



then the internal structure of ( 4.21 ) leads to the Jacobi transformations for an n-degrees of freedom 
systems. It is then a pure coordinate change of variables. As before free motion equations corresponding 
to the undeformed matrix and motion with potential corresponding to modified matrix separate in same 
system of coordinate. 



4.3 Momentum dependent change of variables 

Let us consider functions Bj(X) and Aj(A), which are algebraic functions of entries of the matrix L(X) 
and which satisfy suitable Poisson brackets, for example 

{A j (X),A k (n)} = {B j (X),B k (n)} = 0, 
{^j(A),-Bfe(//)} = Sjk ■ g(X, /i) • (-Bj(A) - B k (ji)) . 

We shall require that zeroes of the functions Bj(X) and values of functions A, (A) at the these zeroes 
define a new system of coordinates. Because functions Bj(X) and Aj(X) are meromorphic functions 
it will also be interesting to consider the corresponding residues. As an example, we present such 
transformation related to Henon-Heiles system [ |32"| . 

Let us restrict to the XXX case. In same way as for the Jacobi transformations we consider the 
block matrices L(A) ( |4.19D with the simplest blocks ( |4.2C| ) Q. 



Lemma 5 The map (u,v) 



-> {x,p) given by 
dk(X) - d k +i(X) 



Pk 



b k (X) - b k+1 (X) 
afe(A) - a k +i(X) 



a G R , fc = 1, . 



A=0 



b k {X)-b k+ i(X) 



fe=i 



A=0 



A=0 



71—1 a 

Eap 

k=l 



(4.22) 



/3 G R , 



1 

^ fe=i 



n-l 
a/3 > y — 



A=0 



k=l 



1 + ^ 



where a k (X) and bk(X) are the entries of matrices Lmn(X, v k , u k ); k = 1, . 
t/ie corresponding determinants, is a canonical transformation. 



n ( j.2L ) and d k {X) are 
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Proof. We fix the variable x\ and the Hamiltonian H — 2fc=i^fe(^ = 0)- 
momentum = a;^} follows from the rs-algebra (2.5) 



Then a corresponding 



ZxkPk = {H,x 2 k } 

d k (X) + d k+1 (X), 



a 



dk(fJ-) - dk+i(p-) 



f( =0.A = 



dk(p) 



■Sfc+l 



2a 



(& fc ( M )-6 fe+1 ( M ))2 

dfe(/*) - 4+i(m) afe(/i) - Ofc + i(/x) 



{(ife(A) + dfc + i(A), 6fc(/x) - 6fc+i(/i)}| 



M =0,A=0 



bk(p-) - b k+1 (fi) b k (n) 
a k (p) - Ofc + i(/i) 



&fc+iO) 



//=0 



6 fc (/x) - &fe+i(/i) 



/x=0 



where we have used the explicit form of the Lmjv matrices (4.20) and the technique developed by 
Sklyanin p8[ . The variable x n is fixed by the condition {x k , x n } — and the corresponding momentum 
p n = {H, x n } is calculated from the rs-algebra and from the Hamilton- Jacobi equations. 
By the proof we have used the following relations 



{b k (X),a j (fx}} 

{b k (X),Cj(ji)} 
{b k (X),bj(^} = 0, 



H — A 
A 



(6 fe (/i)-6 fe (A)), 
(a fe (/x) - a fc (A)) 



0, 



which are determined by the r-matrix only and the one relation {d k (X), dj (/i)} = 0, which is given the 



rs-algebra ( |2.5| ) with s^^-matrices ( j2.35| ). 

We emphasize that it is a non-pure coordinate change of variables and that a free motion on 
manifold and a potential motion on it separate in different systems of coordinates. The canonical 



transformation (4.22) is a technical result as yet, but it is known that the nonlinear equations behind 
the particular systems separated in these variables and in the Jacobi variables are related by a Miura 
map p|. 



4.4 Examples 

We restrict ourselves to the r-matrix algebra of the XXX type and shall consider in detail only the 
Taylor projection ( 2.3C ). 

As an example, we discuss here the Henon-Heiles system and a quartic system of two degrees of 
freedom. Their potentials are polynomials of degree three and four, respectively. The Henon-Heiles 
potential reads (A, B, C are constant): 

H = \{PI+Pl) + \{Ax 2 + By 2 ) + x 2 y + Cy 3 . 

It has been extensively studied both in nonintegrable and integrable regimes. It is known [JFi], ^(J to 
be integrable for the following three sets of parameters 

(i) A = B, C=l/3 

(ii) A, B arbitrary C = 2 (4.23) 

(Hi) 16A = B, C = 16/3 

The quartic potential is described by 

H = \(p 2 x +P 2 y ) + Ax* + Bx 2 y 2 + Cy* . 
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It is known EG] to be integrable for the following three sets of parameters 



(<) 


D 


= 6A, 


C = 


.4 


(«) 


B = 


= 12A, 


c = 


16^4 


(in) 


5 


= 6A, 


c = 


8A 



(4.24) 



Separation of variables of these Hamiltonians takes place in the coordinate systems defined in sections 
4.1-4.3. 



The type (i) systems are related to the Jacobi coordinates and to the block matrices L(X) (4.17 



4.19). Let us consider entries C]\[(X,Uk) in matrices Ljv(A, Uk,Vk) (4.2C) which are the following poly- 
nomials of two variables A and u^: 



c N (\,u k ) 



3=0 



i=0 



N-l N-l-i 

E A * E 44S +1 



(4.25) 



i=0 



We remind that brackets [•] + denote a Tailor projection (2.30). Determinants of matrices 
L^(A, ( 4.20 ) are equal to 

AT 

4 fe) =E/f ) K (4-26) 



A generating function of integrals of motion can be taken as a determinant of the matrix (LIE) 
d(A) = detL(A) = YVk=i ■ Hamiltonians for these systems can be defined by 

d (n-l)N 



H N= dX (n-l)N d W 



(4.27) 



A=0 



and their explicit form up to a constant factor reads 



H 



N 



(k) 
N 



4 



n 
k=l 

T {n) + /3vi n) , (3 e R , 



n N 



k=l 



3 = 1 



(4.28) 



where denotes kinetic energy and vjP is a potential. If higher coefficients /jy of the polynomials 
f( k \\) are the same for all degrees of freedom — fjp for Vfc, j then the hamiltonians (4.27) can 
be rewritten as 



k=l 



A=0 



In the case of the Laurent projection [• ]mn ( |2.29|) Hamiltonians ( 4.28 ) are equal to 



n 



N 



(k) 
N 



fe=l 



(«2- E 

i=-(M-l) 



(4.29) 



(4.30) 



Now by exploiting the internal structure of the matrix L(X) and by using the respective Jacobi 
transformations (4.21) we shall illustrate the general scheme through the cases of two and three degrees 

of freedom. For two degrees of freedom after the transformation u\ = 171+92, u>2 = <7i — 12 the 

(2) 

homogeneous potentials of degrees j = 1, 2, ... N read: 



V, 



(2) _ + 



(2) 
; 

(2) 



= ej qi+ e 2 q 2 ; 

V± 2) +e+ + g|) + ea qi q 2 ; 

5 + 4 (g? + 3«n«2) + [ql + 3gfe) ; 



(4.31) 



V, 



(2) 
JV 



AT J 

EE^Kif + t-^-zrVi 



f(2) \ J-i 
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where = (fffl ± fffi ) and C* are the binomial coefficients. 

For three degrees of freedom and T = J^Pk we can choose the coefficients fffi in such way that 
the Jacobi transformation has the simplest possible form 

Ui = (qi - 2q 2 + q 3 ) , u 2 = -3(g x - q 3 ) , u 3 = (q 1 + q 2 + q 3 ) 
with fffl = 6, fx = 18, fjp = 3. The first homogeneous potentials of degree j — 1,2, ... N are 
Vl 3) = (fP - 3A (2) + A (3) ) gi + (A (3) - 2f^)q 2 + (/« + 3A (2) + fi% 3 ; 

+ (2/f - 4/«)(g l52 + 52 g 3 ) + 2(/ 2 (1) - 9/f + f^)q lQ3 ; 



(4.32) 



y (3) - V {3) 4- V f u,rr> a k a l ■ 

j+k+l=N 



(k) 

where the fjki coefficients are expressed through / and the binomial coefficients CI. For two degrees 



(2) (2) \ ^11 1 

of freedom the Hamiltonians and H\ ( |4.28lj4.31| ) coincides with the Hamiltonians of the Hcnon- 



Heiles system and of the quartic system of type (i). 

Systems of type (Hi) are described by the same block matrices L(\) and by the same function 
/mat (A) as type (i) systems. Systems of type (i) are separated in Jacobi systems of coordinates while 
systems of type (Hi separate after non-pure coordinate canonical transformation of variables described 
in Lemma ^. 

The type (ii) systems are separated in generalized elliptic coordinates. The corresponding matrices 
L(X) were introduced in jl8|, il| a nd have been investigated in detail there. 

Hamiltonians ( 4.23| ) and (4.24) correspond to a special choice of the function /at(A); for the Henon- 



Heiles system 

/ 3 (A) = -i(A 3 -i^A 2 -^ 2 A). (4.33) 

The corresponding particular Lax equations and separation of variables were studied in |l5[ for the 
special choice of the function /at (A), variables of separation for particular quartics potentials were 
discussed in |2^] . Some extensions of these system are known |^0| . They contain more general functions 
./mat (A) and can be investigated within our scheme which describes the infinite family of rational 
potentials. 



5 Conclusions 

We have developed here the rs-matrix scheme for natural finite-dimensional integrable systems con- 
nected to the KdV and the coupled KdV hierarchies. But it applies to other hierarchies of integrable 
nonlinear evolution equations as well. For the hierarchy of the KdV type equations the Lax matrix Mq 



(2.45) is independent of the spectral parameter. The corresponding matrices for the AKNS hierarchy 
and for the SG hierarchy have one pole at infinity or two poles at infinity and at zero correspondingly 
[ p2[ . For applying our scheme in these cases one has to redefine the projector $a: 

= <p{n) ■ Resf =QO + Resfj^ , ip(fi) = fi or ip(fi) = n + - . (5.1) 

M 

These two choices of the function tp(fi) correspond to the AKNS an d to the SG hierarchy if we use 
the second rational parametrization for the r-matrix of XXZ type fl2.20|) for the SG equation. The 
Hamiltonians in terms of the root variables Uk of finite-dimensional systems connected with stationary 
flows for this hierarchies were introduced in J^, |j| . 

In order to apply of our approach to restricted flows of the AKNS hierarchy of equations (2l| we 
have to fix the form of the matrix Mq as: 



M (p) 



[l w 
U —/A 
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and to choose the projector in the form (|5.l| ). 

In the next step one has to consider the undeformed matrix Lq(X) corresponding to this choice of 
Mq and an d to prove that it obeys the r-matrix algebra. Let the matrix Lq(X) has the form 

L (\) = L?(\)+L°(\) 

where we split Lq(X) into the polynomial (pole at infinity) part and the D-divisor part corresponding 
to finite poles. One finds the following sequence of Lg° matrices yielding the prescribed Mo and $A/i 
and satisfying the r-matrix algebra (2.4): 



W) = 

Lg°(A) = 
W) = 



1 




- 


-0 


A 




Pi 


-'0 


A 2 - 




2 


A<?2 


-2 Pl 



{pi,qi} = i 



Mi - 
-x 2 + 



2 P 2 

9192 



{P2,92} = 1 ■ 



The analog of the parametrization related to the Newton representation (4.15) |24j for the recurrent 
construction of these matrices is not known as yet. In this case the Hamiltonians are equal to H = 
^A/i [d(A)], but they do not depend on the spectral parameter since the corresponding residues are 
equal to zero. 

The matrix Mq is more symmetric by its entries in this case and respective integrable systems 



are of non-natural type pl| one can use various solutions of the equations (2.19) for the deformation 
coefficients oij. Systems with such deformations will be investigated in details in the forthcoming 
publications. 

Finally we remark that the rs-matrix scheme presented in this paper applies to the discrete time 
analogs of soliton hierarchies. It yields an rs-matrix description of integrable symplectic maps. 
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